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STRESS WAVES IN TRANSVERSELY ISOTROPIC MEDIA — 


THE HOMOGENEOUS PROBLEM 

Elizabeth R.C. Marques and James H. Williams, Jr. 
Massachusetts Institute of Technology 
Department of Mechanical Engineering 
Cambridge, Massachusetts 02139 


SUMMARY 

The homogeneous problem of stress wave propagation in 
unbounded transversely isotropic media is analyzed. By adopting 
plane wave solutions, the conditions for the existence of the 
solution are established in terms of phase velocities and 
directions of particle displacements. Dispersion relations and 
group velocities are derived from the phase velocity expressions. 
The deviation angles - angles between the normals to the adopted 
plane waves and the actual directions of their propagation - are 
numerically determined for a specific fiberglass epoxy composite. 
A graphical method for the construction of the wave surfaces is 
introduced, using magnitudes of phase velocities and deviation 
angles. The simpler results for the case of isotropic media are 
shown to be contained in the solutions for the transversely 


isotropic media 
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INTRODUCTION 


Most studies on of stress waves in solid elastic materials 
refer to isotropic media. The usual model adopted to describe the 
geometry of such waves is based on the theory of sound waves in 
homogeneous materials as first presented by Lord Rayleigh [l]. 
The same concepts of wave front , dispersion relations and 
directivity functions which are used to describe sound waves are 
also used in the description of stress waves in isotropic solids, 
as can be seen, for instance, in the work of Miller and Pursey 
[ 2 ] and Achenbach [3] • Because an isotropic medium has the 
simplest dynamic behavior among all possible cases of symmetry, 
it can be used as a "reference" when anisotropic cases are 
studied. Concepts (such as phase and group velocities , velocity 
surfaces and wave surfaces ) used to describe the propagation of 
waves in isotropic media have also been applied to anisotropic 
media [4-8] • 

This work is focused on the problem of propagation in 
anisotropic media. The behavior of stress waves in one type of 
anisotropic medium, the transversely isotropic material, is 
described for an infinite medium, establishing the parallel 
between such materials and the isotropic case. Also a graphical 
technique for the construction of the wave surface is introduced, 
giving a better understanding of special geometric features that 
frequently occur during the propagation in filamentary materials. 
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BASIC CONCEPTS 


Definition of Phase Velocity 

The definition of phase velocity is formulated by following 
the path of a small segment of a wave front propagating through a 
medium . The segment is assumed to be sufficiently small that it 
can be approximated by a plane wave segment. The concept is valid 
for any elastic medium , isotropic or anisotropic* 

The phase velocity is defined for periodic or nonperiodic 
waves as [9] 


v n " d / * ( 1 ) 

where d is the distance travelled by the plane segment in the 
direction of the normal n to the wave front in time t (see Figs. 
I.a and l.b ) . The subscript n for the velocity indicates that 
the phase velocity represents the speed of propagation of the 
front segment in the direction of the normal n. 

For periodic waves, the phase velocity can be represented 
in terms of either wavelength and period or frequency and wave 
number as [8] 


v = X/T = co/k (2) 

n 

where X is the wavelength , T is the period of the wave, co is 
the radian frequency and k is the wave number magnitude. 

It is assumed that, in general, the plane wave segment 
travels along a direction that is different from the the wave 
front normal n. If this direction is called r (as for ray of 
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propagation ),the angle between £ and n ia defined as the 
deviation angle [10]. It is known that the deviation angle for 
isotropic materials is zero [9] • 

Also, it can be shown that stress waves in infinite elastic 
media are nondispersive [4*6] . As a consequence , the phase 
velocities and deviation angles are independent of frequency for 
any symmetry configuration of properties of the medium , 
including the isotropic and transversely isotropic media 
considered here . 
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Definition of Group Velocity 

The group velocity is defined as the velocity of propagation 
of energy in a medium [4]» 

For a periodic wave propagating in a medium, the components 

v , of the group velocity in a cartesian coordinate system x, 

Si 

y and z are [9] 

v gi = d cu /dk (3) 

where the k^ 's (i=x,y,z) are the components of the wave 
number vector along the coordinate axes and oj represents the 
radian frequency of oscillation of the wave . The frequency is 
usually represented as a function of the wave number components 
and the elastic properties of the medium ; this representation 
is known as the dispersion relation [ll]. If the dispersion 

relation is known, the derivatives in eqn.(3) can be easily 

calculated. These relations will be derived later from the 
expressions of phase velocity for the media of interest, 
including isotropic and transversely isotropic materials. 

In unbounded elastic isotropic media the elastic properties 
are the same in all directions. For this reason it can be shown 

that the magnitudes of the phase and group velocities are equal 

and the directions £ and £ are coincident [7] (see Fig.I.a). 
These results are to be expected since the system is 
nondispersive . 

For a transversely isotropic medium, the magnitudes of phase 
velocity and group velocity are distinct, and it can be shown 
that the deviation angle , between n and £,has the value [ 1 0 ] 
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( 4 ) 


COS /\ = V /v 

n g 

2 2 21 /2 a 

where v = (v +v +v ) ' . The deviation angle A is 

§ g x sy 

indicated in Fig. 1b * 

The physical meaning of a nonzero deviation angle is that 
energy propagates obliquely with respect to the wave front 
normal. The fact that nondispersive systems can have distinct 
values of phase and group velocities was first emphasized by 
Lighthill [12]. 
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PLANE WAVE SOLUTION OF THE HOMOGENEOUS PROBLEM 


The solution of the equations of motion is derived for the 
specific problem of the propagation of stress waves in 
transversely isotropic media. All the assumptions and results 
also apply to isotropic media ,as they are a subclass of 
transversely isotropic media . In the following sections all 
results are first derived for a transversely isotropic medium and 
then simplified to an isotropic medium. 


Equations of Motion for Transversely Isotropic Medium 

The following assumptions are made for the derivation of the 
equations of motion: 

(1) The medium is homogeneous and has constant density. 

(2) The medium obeys Hooke's law. 

The equations of motion can be obtained for any elastic 
medium by the force-dynamic requirements of a volumetric element. 
Using a cartesian coordinate system, oxyz, as the reference 
system, these equations can be written as [9] 


T + T +T 

xx, x xy,y xz,z 

T + T +T 

xy, x yy,y yz,z 

T + X + T 

xz, x yz,y zz,z 


P u 


,tt 


P V 


,tt 


p w 


, tt 


(5) 


where 


x 

rs 


u,v and 


(r,s=x,y,z) are the normal (r=s) and shear 
(r*s) stresses with respect to system oxyz; 
w are the displacement components of a point in 
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the medium along the directions x,y and 
z, respectively; and 
P is the density. 

The indexes following commas denote derivatives. 

Now, assume a general orthotropic medium in which the 
principal axes coincide with the reference system oxyz. The 
generalized Hooke's law for this medium is [13] 


xx“ 

C 11 u ,x + 

c iz r ,j* 

C 13 w »z 

yy = 

C 12 u ,x t 

c 22 y ,y* 

°23 w »z 

zz~ 

C 13 u » x + 

C 23 v .y* 

C„w 
33 , z 

xz - 

C 44 ( V 

*",x ) 


yz = 

C 55 (v ,x 



xy - 

C S6< u ,y 



the 

°1J are 

the nine 

independent 


( 6 ) 


stifness matrix. 

Here the interest is focused on transversely isotropic 
media, which comprise a subclass of the orthotropic media 
described by eqn.(6), where the number of independent constants 


is reduced to five by the following constraints: 


C 11 =C 22» 


C 13 =C 23’ 


C 44 =C 55 and C 66 - ^ C 11” C 12^ 2 


(7) 


In accordance with eqn.(7), the xy plane is taken to be the 
isotropic plane for elastic properties. 
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The equations of motion, eqn.( 5 )* can be written in terms of 
the displacements u, v and w , using eqns.(6) and (7) for the 
transversely isotropic medium ,as 


C ^u +C._v +C„,w +C t ^(u +v )+C..(u +w ) s p u 

11 ,xx 12 ,yx 13 ,zx 66' ,yy ,xy' 44 » zz ,xz' w ,tt 

C^z-Cu +v )+C„_u +C...V +C.-W +C..fv +w )= Pv .. (8) 

66' ,yx ,xx 12 , xy 11 ,yy 13 ,zy 44' ,zz ,yz' ,tt v J 


C, ,(u +w ) + C..(v +w )+C.,(u +v )+C,,w = p w . , 

44' ,zx ,xx' 44 >zy ,yy 7 13 ,xz yz 33 ,zz ,tt 


Conditions for Existence of Plane Wave Solution 


Now, assume a plane wave solution of the form [ 10 ] 


(u,v,w)= (P ,P ,P ) exp[ 2 JJ_i(xn +yn +zn -v t)l 
' ' ' x y z — J— x y z n J 


where P^, P^, P z are the amplitude components of a particle 

displacement vector along the coordinate axes x, y and z, 

respectively, corresponding to a plane wave with unit normal n; 

n , n and n are the direction cosines of the unit 
x y z 

normal n. along the coordinate axes x, y and z, respectively; v r 
is the phase velocity defined for the direction ri , and X is 
the associated wavelength with the selected n and v r . 

If the assumed solution , eqn.(9), is substituted into the 

t 

equations of motion, eqn.(8), the following expressions are 
obtained [ 10] : 


[ C . n Z +C ce .n 2 +C . .n 2 - Pv 2 ]P +(C. 0 +Cz-^)n n P +(C.^C A A )n n P =0 
L 11 x 66 y 44 z n J x ' 12 66 ' x y y 13 44 x z z 
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< C 12* C 66 )n xVx + t C 66 n / +C 11 n J r 2+C 44 n z 2 - P V n 2] V (C 13 t °44 ) Vx=° 
(C 13 +C 44 )n x n x P x +(C 13 C 44 ) Vz P y* [C 44 ( "x 2+n y 2)iC 33 n z 2 - 


( 10 ) 


The condition for the existence of the plane wave solution 
can be expressed by setting the determinant of the matrix of the 
coefficients of P , P and P in eqn..(lO) equal to zero [10,14] s 

X V 2 


[C„.n 2 +C c ,n 2 
1 11 x 66 y 

+C n 2 -pv 2 ] 
44 z H n J 

(C._+C rr )n n 
'12 66 x y 

^ C 13 +C 44^ n x n z 


(C 4 ,+C ic )n n 
' 12 66 x y 

[C — n 2 + C.„n 2 + C n 2 - pv 2 ] 
66 x 11 y 44 z n J 

( C 13 +C 44^ n y n z 

=0 

( C ,+C . . )n n 
13 44 x z 

(C._ + C . . )n n 
13 44 y z 

[C 44 (n x 2+n 

2 )+C n 2 -pv 2 ] 
Y 33 z K n J 



(ii ) 


Eqn.(ll) is known as Christoff el ' s equation [5] and is 

written here specifically for transversely isotropic materials, 

considering the symmetry conditions and relations between the 

elastic constants of the stiffness matrix given in eqn.(7). The 

expression for the determinant is a cubic equation in v , and its 

solution gives three possible values of v for each selected set 
® * n 

of n , n and n . ( Recall that each unit normal vector n 

x y z 

in the medium is uniquely associated with a set of components n^, 

n and n .) Moreover, the elastic properties are symmetric 
y z 

with respect to the z axis, so the direction cosines for the 
normal ri can be expressed in terms of the angle between n and 
the z axis, as [ 14 ] 
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( 12 ) 


2 2 2 
+ n y = sin 0 

2 2 ft 
n = cos 9 
z 

The advantage of expressing the direction cosines in terms 
of the angle 0 is that the expressions for the velocities can be 
written in terms of a single variable. 


Expressions for Phase Velocities 

By combining eqns.(ll) and (12), the expressions for the 
phase velocities associated with the possible plane wave 
solutions can be written as follows: 

-For a transversely isotropic medium 

(v n ) I =[C 66 sin 2 9+C 44 cos 2 0)/ P ] 1/2 

(v n ) IX =[(C 44 +C 11 sin 2 0+C 33 cos 2 9-l/^ )/2P ] 1/2 (13) 

(v n ) III =[(C 44 +C 11 sln20+C 33 COs29+ ^ _ )/2P ]1/2 

where 

£ =[(C 11 -C 44 sin 2 0+(C 44 -C 33 )cos 2 0] 2 +4(C 13+ C 44 ) 2 (sinecose) 2 
-For an isotropic medium 


(r nh m 

[c 66 /p J 1 ' 2 



[C 66 /P J 1 ' 2 

( 14 ) 

( v n^III = 

[c 11 / p ] 1/2 
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Where for the isotropic medium the following constraint relations 
between the elastic constants of the stiffness matrix have been 
applied: 


C 1 1 =C 22 = C 33 * 


C 12 = C 13 =C 23 * 


C 44 =C 55 -C 66 


= ( C 1 1 ” C 1 2 


)/2 


(15) 


The isotropic medium results may be used as a means of 
checking the more general results for the transversely isotropic 
medium. The results in eqn.(l4) are derived from the 
corresponding velocity expressions for the transversely isotropic 
medium. Observe that the values obtained for the phase velocities 
in eqn.( 14 ) are the widely known values of velocities of shear 
waves, ( v n )j and ( v n )jj» an< * longitudinal waves, (v which 

in terms of elastic properties can also be represented by 
(G/P)^ 2 and (E/ P )"*^ 2 , respectively. The values G and E • are 

the shear and extensional elastic moduli , repectively. 


Particle Displacement Vector Components 


By using the computed values of the phase velocities , the 

amplitudes of the particle displacement components associated 

with each of the phase velocity vectors can be found . This can 

be achieved by solving eqn.(10) for the displacement amplitude 

components P^, P and P z for each of the phase velocities. 

Because at each point there are three velocities for each 

direction defined by n, there are also three displacement 

amplitude vectors for each direction n, namely, (p)„=(p .P .p ) 

— J 'I'x’y* z'l’ 

’^y * an< * (£ ) iii = (^ x * ^y » ^ 2 ) in corresponding to 
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(v ) T , (v )__ and (v )___, respectively. The three displacement 
n n ii n j. x x 

vectors are normal to each other but usually neither normal to 
nor coincident with the direction _n for a transversely isotropic 
medium. 

In the homogeneous problem, there is no unique solution for 

the displacement amplitude components. Each phase velocity (v )_, 

n I 

(v n )jj and ( v n )iII ^ as a corresponding displacement vector 

(P)j, (Vij and (P)jjj* respectively , which represents a 

mode of vibration. Any multiple of the three displacement 

vectors (,P)j» (P)jj and (P)jjj is also a solution. 

As there are no initial conditions in the homogeneous 

problem, the components P ,P and P obtained by solving eqn.(lO) 

x y z 

for each (v^)^ can be expressed in terms of an arbitrary 

constant. Thus, the useful information that can be extracted 

from the solution of eqn.(10) for the displacements consists of 
the direction of the particle displacement vector with respect to 
the coordinate axes. The magnitude of the displacement vectors 
(P)^ is expressed in terms of the components as 


(P) 


i 




2 ) 1/2 


i 


1 = 1 , 11,111 ( 16 ) 


Then, the direction cosines of the displacement vectors with 
respect to the xyz coordinate system are [10] 



[n * 2 < C 12* C 66 )]1/2 

B [ pv ^-(n ^+n 2 )C^----n ^C..] 
L H n v x y ' 66 z 44 
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(17) 


P 

_ 2 L 


[\ 2 ' C 12 + C 6 S ) J 1/2 


B C pv n 2 - (n / +n y 2)C 66- n z 2 C 44 ] 


t n z 2 < C 13 +C 44 )2/(C 12* C 66 ” 1 /2 


O ? 9 9 ^ T 

Pv n ' (n x tn y ) C 44- n z C 33* n z (C 13 +C 44^ 

<° 12 * C 66 > 


where 


B = 


(o » 2 * n v 2) (C 12^ C 66 ) 
[ ^ v n 2 - ( "z 2+n /> C 66- n z 2c 44 ]2 


K 2(C 13* C 44 )2/ ' (C 12 tC 66 )] 


) 1 / 2 


Pv 2 -(n 2 +n 2 )C -n 2 C +n 2 (C., + C J 21 2 

n v x y ' 44 z 33 z 13 44 

< C 12 +C 66> 


Observe from eqn.(l7) that the sum of the squares of the 

ratios P /P, P /P and P /P is equal to one, meaning that only 
x y z 

two of the three components can be independently determined. 

The isotropic medium is a particular case of the 
transversely isotropic medium and where the particle displacement 
directions are either normal to or parallel with the direction 
n. In this case, the displacement vectors are aligned with or 
perpendicular to the direction n. 
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VALUES OF ELASTIC CONSTANTS OF THE STIFFNESS MATRIX 


In order to demonstrate the application of the plane wave 
solution for the homogeneous problem of wave propagation in 
infinite media , numerical values are fixed for the elastic 
constants of the stiffness matrix, representing the two media of 
interest: transversely isotropic and isotropic materials. 

Transversely Isotropic Material 

For the transversely isotropic material, there are five 
independent elastic constants in the stiffnes matrix that will be 
taken as C^, C^, C^, C^ and C^. Note that = C^^-ZC^ . 

In this case, the material is a unidirectional fiberglass epoxy 
composite having a resin content of 36 percent by weight and a 
density of 1850 kg/m^. The elastic constants of the stiffness 
matrix were determined experimentally by ultrasonic techniques 
from samples of the material as [15] 

C = 10.581 x 10 9 N/m 2 
1 1 

C 13 = 4.679 x 10 9 N/m 2 
C„ = 40.741 x 10 9 N/m 2 

DJ 

C = 4.422 x 10 9 N/m 2 

44 

G„ = 3.243 X 10 9 N/m 2 

OO 

The constant was calculated from the values of C ^ and 

c above as C..= 4 . 098 x 10 9 N/m 2 . See Fig. 2 for the relation 
66 1 2 

between the coordinate axes and the material's orientation. 
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Isotropic Material 

For an Isotropic material, there are two Independent elastic 

constants In the stiffness matrix that will be taken as C„. and 

1 1 

C^ 2 . For the representative Isotropic material, E glass was 
chosen. 

The values of the corresponding elastic constants are 
C 11= 82.658 x 10 9 N/m 2 and C 12 = 23.31 x 10 9 N/m 2 with a density 
of 2540 kg/m^ [16] . This material in fiber form is one of the 
constituents of the transversely isotropic material described 
above . 
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POLAR DIAGRAMS OF PHASE VELOCITIES 


The values of the phase velocities obtained from eqns.(l3) 
and (14) » using the numerical values of the material properties, 
can be presented in polar diagrams. Such polar diagrams of phase 
velocity are called velocity surfaces [5]. 

The velocity surfaces are shown in Figs. 3 and 4 for the 
isotropic material and the transversely isotropic 
material , respectively. One quarter of each surfaces' intersection 
with the xz plane is presented . The surfaces are designated as 
V(SH), V(SV) and V(P) corresponding to the velocities indexed as 
I, II and III, respectively. 

Among the directions defined by the normal n used to 
calculate the points of the velocity surfaces, several directions 
were selected and presented in. Figs. 3 and 4* These directions 
are represented as the radii from the origin of the coordinate 
system to the surfaces. The particle displacement vectors 
corresponding to these directions are represented by the arrows 
and dots , giving an indication of the relative position between 
the particle displacement direction and the direction of the 
normal to the wave front n. The relative position of the particle 
displacement direction and the corresponding direction of the 
normal ri (radii from the origin ) can be used to identify the 
modes of propagation. If these two directions are coincident, the 
mode is purely longitudinal. If the two directions are 
perpendicular , the mode is purely transverse. 

Observation of the V(P) surface in Fig. 3 shows that 
directions ri and particle displacement directions are coincident 
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for all directions in the medium, indicating the existence of a 
pure longitudinal mode of propagation . The surfaces V(SH) and 
V(SV) in Fig. 3 show that directions n and particle displacement 
directions are normal to each other for all directions in the 
medium, indicating the existence of two pure transverse modes of 
propagation. So, for isotropic materials there are three possible 
pure modes of propagation, one longitudinal and two transverse. 

For the transversely isotropic material, there is one 
possible pure transverse mode propagating in all directions of 
the medium, corresponding to the V(SH) surface as can be seen in 
Fig. 4 . Pure longitudinal modes occur only at the intersections 
of the V(P) surface with the coordinate axes (x, y and z) where 
the particle displacement direction coincides with the direction 
n. For the V(SV) surface there are three directions where pure 
transverse modes occur, namely, the directions of the coordinate 
axes x, y and z. (In Fig . 4 only the intersections with the x and 
z axes are shown) . 



According 

to the 

convention 

adopted 

for the designation 

of 

the 

surfaces , 

Table 

1 shows 

the 

modes of vibration for 

points 

along 

the directions 

of the 

coordinate 

axes and the 

surfaces 

that 

contain 

the indicated modes 

. Other 

designations 

for 

the 

surfaces can 

be found 

in the 

literature 

[10,14,17] , 

as , 

for 

instance, "transverse" 

for the 

SH 

waves , 

"qua si- transverse " 

for 


the SV waves, and "quasi-longitudinal” for the P waves . 
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DISPERSION RELATIONS AND GROUP VELOCITIES 


After the phase velocities for each direction n are 
obtained, the group velocities can be determined by the use of 
the dispersion relation. Combining the definition of phase 
velocity, eqn.(2), and the expressions obtained for the 
velocities, eqns.(l3) and (14) > the dispersion relations can be 
written as follows: 

-For a transversely isotropic medium 

( « ) x = {t C 66 ( k x 2+k y 2)+C 44 k 2 2]/P}1/2 

C to ) IZ = {[C 44 k 2 + C 11 (k/ + k y 2 )+C 33 k z 2 -KT?]/2 P } 1/2 (18) 

( 0) ) 1Z1 = {[C 44 k 2 + C 11 (k x 2 + k y 2 ) + C 33 k z 2 + l/7^]/2 P} 1/2 

-For an isotropic medium 

( a) ) = [C fifi (k 2 + k 2 +k 2 )/P]^ 2 

v I L 66 x y z ' ' J 

( “ )u - •f C 6 6 (lt x 2 * k y 2+lc z 2 )/ e J 1/2 (19) 

( CO ) IIr = tC^U^.k/.k/j/P ] 172 

where the square of the wave number magnitude is represented by 

k 2 = (k 2 +k 2 +k 2 ) (20) 

x y z 

and the relations between the components of jc, namely, 

k 2 + k 2 = k 2 sin 2 0 , k 2 = k 2 cos 2 0 (21) 

x y z 

were used (as in eqn.(l2) for n) . 
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From the dispersion relations in eqns.(l8) and (19) and the 
definition in eqn.(3)> the values of the group velocity 
components along the coordinate axes can be determined . The 
expressions for the group velocity components in the xz plane are 
as follows : 

-For a transversely isotropic medium 

(v )_ = n C c , [ ( C-^-n 2 + C.,n 2 ) P ] -1 / 2 
v gx'I x 66 1 66 x 44 z 

( V’l * V 44 [(C 66 n x 2+C 44 n * 2) p ] ' 1/2 

* v gi^II = ^ n x ^ C 44 +C 1 1 E ^ n z^ C 1l' C 44^* C 1l" C 44^ n x * 

^ C 4 4” C 33^ n z^^ + ^ n x n z^^13 +C 44 )2 ^ 

( 22 ) 

(v gz ) n‘ e " 1/2n z( C 44- C 33 )[(C ir C 44 )n x 2+ 

(C 44- C 33 )n z 2]+2n ^/ (C 13 +C 44 )2t 

(T «x ) IIl' (1/CHn * (0 *4 +C 11 )+ ' ' 1/2n x< C 11- C 44 )t(C 11- C 44 )n x 2t 

^ C 44 _C 33^ n z ^ +2n z n z ^ C 13 tC 44* * 

^' l Sz'lII = * 1/,C ^ n z^ C 44 +C 33^' f S ^ n z^ C 44 _C 33^* C 11 _C 44^ n x + 

(C 44- C 33 )n x 2]+2n 2 n x 2(C 13 +C 44 )2) 

where B={ t C^Cn/en^ 2 )* c i 1 ***<>„* *-]fT 1/2 P ) 1/2 
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C * (tC 44 (n x 2+n / ,+C 11 n /* C 33 n / t l /T 1/20 )1/2 


-For an isotropic medium 


2._ 2, n -l/2 


^gx^ ( V ) II = n x°66 [C 66 p (n x +n z ^ 

( v gzh = ^gz^II = n z C 66 ^ C 66 P (n x 2+n z 2 ^ 1/2 

(v ) TTT = n C. . [C,„ P (n 2 +n 2 )] -1 / 2 
v gx III X 11 L 11 ' X z /J 

(v ) TTT * n C . A [C,* P (n 2 +n 2 )]" 1/2 
v gZ III Z 11 1 11 X z /J 


( 23 ) 


Observe that the magnitudes of the group velocities can be 
obtained from the corresponding components in the xz plane as 


(v ) = [(v 2 +v 2 ) 1 ^ 2 ] 


(24) 


For an isotropic medium, the magnitudes of the group 
velocities are identical to those of the phase velocities as 
defined in eqns.(l4) . 
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CALCULATION OF DEVIATION ANGLES 


The purpose of calculating the deviation angles is to 
determine the actual direction r of propagation of plane wave 
segments having normal n. The deviation angle can be expressed 
in terms of the phase and group velocity magnitudes as given in 
eqn.(4) . 

Because z is the symmetry axis in the transversely 
isotropic medium under consideration , the directions of phase 
velocity and group velocity (n and r, respectively) are expressed 
in terms of their angles with respect to this axis . 

For the phase velocity direction , the angle 0 with 


respect 

to the z axis ( 

see Fig. 1b ) is 

defined 

by the direction 

cosines 

of n as given in 

eqn.(l2) as 





< tan 9 >1, II, HI* [(n x 2+n y 2)1/2 

/ n z ] I, 

II, III 

(25) 

The 

angle 0' of 

the direction r 

( for 

which the 

group 

velocity 

is defined) 

with respect 

to the 

z axis 

can be 


expressed in terms of the group velocity components as 



The indexes I, II and III in eqns.(25) and (26) denote that 
the same equations apply to each of the velocity sets. Moreover, 
due to symmetry, n and v may be taken as zero ,and the 

y & y 

calculation can be done for the xz plane . 

The general expression for the deviation angle is then 


A = 9 


(27) 
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The specific algebraic expressions, expressed in terms of 
elastic constants for the deviation angles of transversely 
isotropic media, are not reproduced here. For isotropic media 
the deviation angles are zero. It is interesting to note that 
algebraic manipulations show that the deviation angle is 
independent of the frequency. And, since the phase velocity is 
also independent of frequency (see eqns.(l3) and (14)), so is 
the group velocity, which is another way of expressing the 
nondispersive character of the elastic media under study. 

The calculated deviation angles for the transversely 
isotropic medium under consideration are shown in Figs. 5 a , 5 b 
and 5c , for the SH,SV and P modes, respectively. 
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GRAPHIC CONSTRUCTION OF WAVE SURFACES 


In this section the graphic construction of the wave front 
will be described, focusing on the transversely isotropic 
material. The quantities needed for this construction are the 
magnitudes of the phase velocities ,the directions _n for which 
these velocities are defined and the deviation angles. 

In other words , the locus of the points of equal phase in 
geometric space for unit time will be determined, supposing that 
plane waves passed through the origin at time t=0 in all possible 
directions defined by their normals n. The construction is 
described for the xz plane but it is valid for any plane 

containing the symmetry axis z. 

The following steps must be followed for the construction (see 
Fig.(6))s 

(1) Choose a direction n (n , n ) 

X z 

(2) Draw a line from the origin of the coordinate system 
in the direction ji and having a scaled length equal to 
v fl . Since time is taken as unity, the segment represents a 
distance . 

(3) From the line with direction n, mark the deviation 
angle , if positive below ,if negative above the line 
having the direction n. The direction so determined is £ 
which is the direction of the the group velocity v^ . 

( 4 ) From the tip of the line of length v r , draw a line 
perpendicular to it. 

(5) The intersection of the perpendicular line with the 
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direction r determines a point P on the wave surface. 

Thus, point P represents the location at time t=1 sec of a 

point of the wave front which at time t=0 passed through 

the origin with velocity v in the direction of n. 

n 

If the construction described is repeated for each set of 
calculated phase velocities assuming all possible directions _n 
of the wave normal in the medium , the wave surfaces are 
obtained. The construction steps of the wave surfaces W(SV) and 
W(P) corresponding to velocity surfaces V(SV) and V(P) 
respectively, are shown in Figs. 7 and 8 for some directions in 
the representative transversely isotropic medium. In these 
figures, the primed points are points of the extremities of 
the segments of length v n as described in step 2). The umprimed 
points are the corresponding points on the wave surface. 

The construction of the wave surface W(SH) corresponding to 
velocities V(SH) is not shown since the deviation angles are 
small and the final shape of W(SH) is very similar to the shape 
of V(SH). The three wave surfaces obtained for the 

representative transversely isotropic medium are shown together 
(positive x-z quadrant only) in Fig. 9 • 

For the isotropic medium no construction is made since the 
deviation angles are zero. For such a case, the wave surfaces are 
identical to the corresponding velocity surfaces. 

It is important to emphasize that the wave surfaces are the 
actual geometric positions of wave fronts for the three possible 
modes of propagation. Because plane wave solutions are assumed to 
generate these surfaces , the wave surfaces represent the envelope 
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of all possible plane waves emanating from the origin of the 
coordinate system in all possible directions, as if such a 
source existed at the origin. 
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SPECIAL FEATURES OBSERVED IN THE WAVE SURFACE W(SV) 

As it can be observed in Fig. 7, for the representative 
transversely isotropic medium there is a folding of the W(SV) 
wave surface due to an inversion of the sign of the deviation 
angle 

To explain the geometric shape of the wave surface W(SV), 
assume that the wave front is constituted of many small segments 
of plane waves having the same length and passing through the 
origin at time t=0. Assume that for an infinitesimal interval of 
time these segments of plane wave are still very close to the 
origin and form a polygonal line that can be inscribed in an arc 
of circumference. Then immediately thereafter, the segments start 
to deviate from their normals, following the directions r defined 
by the deviation angles .If it is assumed that the segments' 
extremities remain connected physically , the wave front must be 
continuous in time as these segments will be stretched as they 
travel away from the origin. 

The deviation angles, Fig. 5b, can be associated with a 
degree of stretching of the original straight segments . The 
larger the deviation angle , the larger the degree of stretching 
undergone by the segment. Recall that the angle 0 represents 
the angle of the normal n (to each segment) with respect to the z 
axis. 

Observe that the curve of has two maxima, one 

positive at 0 = 33° and one negative at 9 = 77.5°, and passes 
through zero at 0 = 60°. If the curve is divided in four 
regions as shown in Fig. 10b ,it can be observed that in region 1 
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the segments are stretched increasingly until the first maximum 
at 0 = 33°. Then for region 2 the segments are less and less 
stretched until the deviation angle reaches the value 
zero, represented by point C in Fig. 5. a. The segment 
corresponding to the zero deviation is not stretched. For region 
3, again the segments are stretched from zero to the second 
maximum at 0 = 77*5° and finally in region 4 they are less and 
and less stretched until reaches zero for the segment whose 
normal is located at 90° with respect to the z axis. 

In the model of the polygonal line constituted of straight 
segments, if points A and E, Fig. 10. a, are the intersections of 
this line with the coordinate axes, for an infinitesimal time 
after zero time , these points remain on the coordinate axes since 
the corresponding deviation angles are zero. Later ,at a time 
t equals unity, these points will be at positions A' and E', 
respectively. Considering that at the same time all segments are 
supposed to remain connected ,the segments originally in regions 
1 and 2 (Fig. 10. b) will undergo an effect of pulling away from 
the z axis; the segments originally in regions 3 and 4 on the 
other hand , will undergo an effect of pulling away from the x 
axis, moving closer toward the z axis. 

It has been determined that the segments which are the tips 
of the cuspidal edges B',D' of the wave surface W(SV) have 
normals n at angles of 54*02° and 75*06°, respectively, with 
respect to the z axis. For time t equals unity, the location of 
these two segments is defined by the lines passing through the 
origin and the points B' and D', respectively, at 61.86° and 
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41.93° with respect to the z axis (see Fig. 10a ). 

The originally equal segments that at time t equals unity 
are located inside the region limited by the lines OB’ and 0D' can 
be associated with the amount of energy that is concentrated in 
this region. If it is assumed that each segment that passes 
through the origin at time t=0 carries the same amount of energy, 
the percentage of the total number of the segments contained in 
the region of interest at time t equals unity is also the 
percentage of energy contained in the region at any other time 
t. Using this concept, it was found numerically that 64 . 4 $ of 
the total energy propagates through the region located between 
the lines OB' and 0D • where the included angle B'OD' is 20.93°. 
This calculation was made for a total number of 180 segments. 

Concerning the geometry of propagation of the entire wave 
front, it can be said that if an observer stands at a position 
with respect to the coordinate axes inside the region defined by 
lines OB' and 0D ' (for instance, at position P in Fig. 10. a ), 
such an observer will see three portions of the same front W(SV) 
passing through him at different times . 

It must be stressed that the special feature described above 
is valid for the fiberglass epoxy composite used in this 
numerical example. The existence of this special behavior is not 
a common characteristic of all transversely isotropic media and 
depends on the numerical values of the elastic constants , so 
such features must be analyzed and determined for each specific 
material . 
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CONCLUSIONS 


The problem of the geometry of propagation of waves in 
transversely isotropic media was studied. The solution of the 
homogeneous problem of motion for a fiberglass epoxy composite 
was analyzed and the numerical values of the resulta were 
presented . 

The solution of the homogeneous problem of motion in an 
infinite medium provided a description of the behavior of plane 
stress waves in terms of velocities and paths of propagation 
Using these concepts, a graphical method for the construction of 
the wave fronts was introduced . 

As described above, the following conclusions can be stated: 

( 1 ) The graphical technique is found useful in the construction 
and the physical interpretation of the wave surfaces . The 
technique was applied to transversely isotropic materials but 
the concept can be extended to other cases having anisotropic 
symmetry . 

(2) The wave surfaces corresponding to a fiberglass epoxy 
composite material were constructed .Special features were found 
in the W(SV) surface. 

(3) It was found that energy carried by the SV waves in the 
fiberglass epoxy composite material is concentrated along certain 
preferential directions. This phenomenon can be explained in 
terms of the spreading geometry of plane wave segments . 

The conclusions above suggest some possible applications and 
recommendations for ultrasonic NDE of transversely isotropic 
materials . For instance, the special features of the wave 
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surface W(SV) for the fiberglass epoxy material can be used as 
an auxiliary criterion for the placement of transducers when 
experiments are designed . For example, the locations where the 
energy is focused (angle B'OD' in Fig. 10 ) can be selected as the 
receiving region since the signals within this region are 
expected to be stronger due to larger displacement amplitudes. 
Also, certain locations in the medium can be selected as 
receiving points where up to five different wave front arrivals 
may occur, when the three modes of propagation are considered. 

The knowledge of the geometry of the spreading can also be 
used for crack detection since the number of expected arrivals 
can be computed for each direction in the medium. Thus, the 
absence of the arrival of one of the expected portions of a wave 
front might indicate that an obstacle was in the path of that 
portion of the front. Such an obstacle might indicate the 
presence of a flaw. 
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Table 1. Designation of Velocity Surfaces and Modes of Vibration 
at Intersections with Coordinate Axes. 
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g.2 Position of material axes with respect to coordinate 
axes for transversely isotropic medium. 




Component, of phase velocity in z direction (m/s) 



Fig. 3 


Phase velocity polar diagram for isotropic E glass 
fiber material for positive x-z quadrant. (Arrows 
represent directions of particle displacements.) 




Component of phase velocity in z direction (m/s) 





Deviation angle Deviation angle Deviation angle 



a) Deviation angles 
for V(SH) 



Deviation angles 
for V(SV) 



c) Deviation angles 
for V(P) 


Fig . 5 Deviation angles ( ) T TT TTT between phase 

velocity and group 1,11,111 ve iocity directions 

for phase velocities V(SH),V(SV) and V(P). 
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Fig. 7' Construction of wave surface W(SV) from 
corresponding velocity surface V(SV) and deviation 
angles (A ) . ^ for fiberglass epoxy composite , 
for positive x-z quadrant. (Points 1 through 13 
are on wave surface : points l 1 through 13’ are on 

velocity surface.) 
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Construction of wave surface W(P) from 
corresponding velocity surface V(P) and deviation 
angles (/\J _ _ for fiberglass epoxy composite , 
for positive 1 x-z quadrant. (Points 1 through 10 
are on wave surface : points l'through 10' are on 

velocity surface.) 


Distance traveled by wavefront in z direction 
in I second, z(m) 





a) Wavefront assumed shapes for time t = 0 ( ABCDE ) 
and t=unity (A'B'C'D'E') . 



0 20 40 SO 90 


Angle of n with respect to z axis, 0 (°) 

b), Deviation angles (A) showing regions and 
characteristic points! 


Fig.lo Model of geometrical spreading for SV wave. 
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